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COMMON APPROACHES IN DESCRIPTION 
OF ORDINARY LIQUIDS AND HADRONIC MATTER 


This work is an attempt to give a brief overview of the implementation of the statistical ther¬ 
modynamics to hadronic matter. The possibility to use the hydrodynamic approach for develop¬ 
ing the physical model of the formation of exotic structures in the head-on intermediate-energy 
heavy ion collisions is discussed. That approach is shown to be able to provide simple analytical 
expressions describing each step of the collision process. This allows for extracting the data 
concerning nuclear matter properties (surface tension, compressibility, etc.) from the proper¬ 
ties of observed fragments. The advantages of the thermodynamic analysis of phase trajectories 
of the system in heavy ion collisions are discussed. Within the thermodynamic approach, the 
method to evaluate the curvature correction to the surface tension from the nuclear matter 
equation of state is described. The possibility to use statistical thermodynamics in the studies 
of hadronic matter and quantum liquids is discussed. 

Keywords: heavy ion collisions, equation of state, hydrodynamic instability, phase diagram, 
Tolman length, thermodynamics, Skyrme parametrization. 


1. Introduction 

The study of the properties of hadronic matter and 
the phase transitions in it has been the subject of 
many studies both theoretical and experimental in 
recent decades. Such an interest is explained both by 
the applied importance of results for the energy pro¬ 
duction industry, radiation medicine and by their fun¬ 
damental importance, because the many-body sys¬ 
tems, in which the excitation energy is thermally dis¬ 
tributed among many particles and in which there ex- 
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ist long-range Coulomb and short-range forces, up to 
now have no satisfactory description. The description 
of nuclear matter (NM) within statistical thermody¬ 
namics adopted in the recent decades is very similar 
to that for ordinary liquids The physical rea¬ 

son for that is the similarity between the molecular 
forces and nuclear-nuclear interactions in respect to 
their dependence on the distance. Notwithstanding 
the tremendous difference in the energy and space 
scales, the nuclear and van der Waals forces represent 
a very similar behavior. In both cases, the attraction 
between particles is replaced by their repulsion at a 
small interaction range. J.D. van der Waals suggested 
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his famous equation in 1873. A hundred years later, 
his finding was fruitfully used to describe the proper¬ 
ties of the nuclear matter unknown to the nineteenth 
century scientists. It is quite natural to expect that, 
when investigating both ordinary liquids and nuclear 
matter, it is reasonable to use similar fundamental 
concepts of the modern thermodynamics, statistical 
physics, phase transitions theory, and the liquid state 
theory accounting for the quantum nature of nuclear 
matter. 

Among all the studies of nuclear systems, of high 
importance are those devoted to the processes ob¬ 
served in heavy ion collisions (HIC). Such studies are 
valuable from the fundamental point of view, as it 
is the only available experimental method (under the 
Earth conditions) to obtain data about the nuclear 
matter at extreme temperatures and densities. As for 
the applied importance, the results are needed to pro¬ 
ceed in developing the new technologies of the nu¬ 
clear waste utilization, which is an up-to-date task 
accounting for the rapid development of the nuclear 
energy production industry. The other important ap¬ 
plications can be found in medical physics and bio¬ 
physics. When the beams of high-energy ions pass 
through the matter, they leave the trace of the frag¬ 
ments produced in the collision, causing the changes 
in structures and thermodynamic properties of bio¬ 
logical objects and liquids. Therefore, the clear un¬ 
derstanding of the fragments formation mechanisms 
can help in biophysical studies and in developing the 
new technologies in the medical physics. 

Nowadays, it could be seen that, in the field of 
studies devoted to the HIC, the progress in experi¬ 
mental researches is much faster than in theoretical 
ones. Unfortunately, there are yet no complete the¬ 
ories of such important phenomenon as the nuclear 
multifragmentation, no satisfactory equation of state 
(EOS) for hadronic matter, and no clear understand¬ 
ing of phase transitions in nuclear matter. Therefore, 
the progress in the theory able to describe the basic 
thermodynamic and structural properties of hadronic 
matter, phase transitions in it, and behavior of hot 
nuclei with high excitation energies stands among the 
most important problems in modern physics. 

To recall the stage, we would like to present a brief 
overview of the problems that exist in hadronic mat¬ 
ter physics and can be studied with the help of statis¬ 
tical thermodynamics and the theory of phase tran¬ 
sitions well developed for ordinary liquids. 


1.1. Nuclear matter equation of state 

The knowledge of the EOS of nuclear matter is one of 
the fundamental goals in nuclear physics, which has 
not yet been achieved It should describe the 

fundamental properties of NM 


P = P{p,T,6), 


( 1 ) 


where P is the pressure, p is the density, and S is 
an asymmetry parameter, at any point of the phase 
diagram of NM (Fig. [T]). The EOS of NM is an im¬ 
portant ingredient in studying the properties of nu¬ 
clei at and far from stability, studying the structures 
and the evolution of compact astrophysical objects 
such as neutron stars and core-collapse supernovae, 
applied nuclear physics, etc. The saturation point of 
EOS, i.e., the energy as a function of the matter den¬ 
sity, for the symmetric NM (SNM) at zero tempera¬ 
ture (T = 0), is well determined from the ground-state 
properties of nuclei, such as binding energies and cen¬ 
tral matter densities, by extrapolation to infinite NM 
ll|. At the same time, one can find the vast majority 
of different parametrizations and different theoretical 
models that were built in attempt to describe the NM 
and finite nuclei in a wide range of the external pa¬ 
rameters. All of them were built for the case of special 
assumptions and are lacking a predictive power (l^ . 

The things go even worth, when one changes to the 
asymmetric NM and the pure neutron matter. Exist¬ 
ing models give different results for the symmetry en¬ 
ergy of asymmetric NM [1^ , its slope, and magnitude 



Fig. 1. Phase diagram of NM E3 


ISSN 2071-0186. Ukr. J. Phys. 2015. Vol. 60, No. 8 


709 










K.V. Cherevko, L.L. Jenkovszky, V.M. Sysoev et al. 



Fig, 2. Left window: Density dependence of the nuclear symmetry energy Esym{p) from SHF with 21 sets of Skyrme interaction 
parameters fisl . Right window: Same as left panel from the RMF model for the parameter sets NLl, NL2, NL3, NL-SH, TMl, 
PKl, FSU-Gold, HA, NLp, and NLp5 in the nonlinear RMF model (solid curves); TW99, DD-MEl, DD-ME2, PKDD, DD, 
DD-F, and DDRH-corr in the density-dependent RMF model (dashed curves); and PC-Fl, PC-F2, PC-F3, PC-F4, PC-LA, and 
FKVW in the point-coupling RMF model (dotted curves) [l^l 



n (fm^l n |fm^l n |fm^| 


Fig. 3. Energy per nucleon in SNM -E'o(n) (left panel), the 
symmetry energy Esin) (middle panel), and the energy per 
nucleon in a nonsymmetric NM (NNM) (B-equilibrated and 
charge neutral) for different models (right panel) [H 

at the saturation densit y a nd, the experimental con¬ 
straints are quite wide [1^. It should be noted that 
the symmetry energy behavior is very different for 
different EOS parametrizations in the high- and 
low-density regions (Fig. [^I- 
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The similar uncertainty is observed for a symmetric 
NM (SNM) in the high-density region (Fig. [3]). 

Nowadays, extracting the information on the nuc¬ 
lear EOS, in particular at high baryon densities, 
is restricted to observations of astrophysical com¬ 
pact objects and studies of hot nuclear systems cre¬ 
ated in high-energy proton-induced reactions or in 
intermediate and relativistic energy heavy-ion colli¬ 
sions (HIC). 

Therefore, HIC can shed light on what are nuclei 
made of and explain the properties of the matter that 
form the nuclei. Such experiments can be treated as a 
kind of materials science similar to the P-V-T studies 
in the physics of ordinary liquids. Hence, it is quite 
natural to apply the well-developed theories of ordi¬ 
nary liquids to the description of HIC. Understanding 
the behavior of NM in the HIC is of great relevance 
either for theoretical nuclear physics or for nuclear as¬ 
trophysics, when studying the supernova explosions 
and the properties of neutron stars (l8l - Eol |. 

2. Heavy Ion Collisions 

The idea of HIC experiments is to collide two nucleus 
or light particles with nuclei at a high enough energy 
to break them into pieces. Introducing then an ap- 
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propriate model of the process involved, it should be 
possible to extract some data on NM properties and 
its behavior at high densities and temperatures. The 
basic mechanisms of the formation of fragments in 
the HIC experiments are the following 

• spallation, when only one heavy fragment with 
its mass close to the target mass is formed (to = 1); 

• fission that goes usually with the two heavy frag¬ 
ments produced with masses in the interval around 
half the target mass A ^ ^; 

• multifragmentation that is the process which 
leads to the formation of several fragments with in¬ 
termediate masses. Usually with A < 50. 

The last one is especially important in studies of 
the NM EOS. Disintegration, when a bigger nucleus 
breaks into one or several nuclei and some nucleons, 
is one of the two basic mechanisms in the Nature, 
by which a nuclei can be formed. The appropriate 
way to study this phenomenon is the experiment with 
proton-nucleus and nucleus-nucleus collisions. When 
such a collision is observed at a high enough en¬ 
ergy (higher than the threshold value Eth ~ 2- 
4 MeV/nucleon (IJ), there is a high probability 
of the production of intermediate-mass fragments 
. The nature of the production of intermediate- 
mass fragments in such experiments is of great in¬ 
terest for the consequent understanding of the NM 
properties and for the progress in the studies of nu- 
24-^. Thus, HIC have been intensively 


clear EOS 


Initial State 



Fig. 4- Schematic picture of a fragmentation reaction, in 
which a given initial channel may develop into many different 
fragmentation channels during the dynamical evolution IstII 

important challenge is that the heavy ion collision ex¬ 
periments are inclusive and will remain so in the fore¬ 
seeable future. That means missing some information 
about the physical processes going on in the system. 
Only the degree of inclusiveness may change if more 
sophisticated counter arrays are introduced 0. 

During the last decades, a wealth of data have 
been accumulated about the multifragmentation phe¬ 
nomenon in proton-nucleus and nucleus-nucleus col¬ 
lisions 

m 


studied in the last years 25, 3^ . It should be noted 
that, in a composite system formed in the collision 
process from the projectile and target nuclei, nuclear 
matter can be highly compressed and be at high tem¬ 
peratures during the early stage of the reaction. At 
the time of the formation of intermediate-mass frag¬ 
ments with 3 < Z < 20, their characteristic prop¬ 
erties such as the excitation energy and isotopic dis¬ 
tributions are governed by the characteristics of the 
break-up source such as the temperature, density, and 
N/Z ratio. Therefore, intermediate-mass fragments 
may provide a unique probe to study the reaction 
mechanism and hot nuclear matter properties 0 - 
At the same time, in order to have the complete pic¬ 
ture of the situation in a field, one should bear in mind 
that the colliding system is over a large time span of 
the reaction out of global and even local equilibrium 
0. That means the presence of the non-equilibrium 
effects all over the compression phase, where one es¬ 
sentially intends to study the NM EOS. The other 


T3,0[3i, [ 3 ^. Different models (statistical 
percolation |35l |i and different mechanisms 
(liquid-gas phase transitions, spinodal decomposition, 
sequential evaporation from the expanding-emitting 
source, etc. [36|) are being used for explaining the ex¬ 
isting data. As for the initial highly non-equilibrium 
stage, the most common approach is using the trans¬ 
port models for its adequate description. One of the 
important points in that case is allowing the model 
for the possibility for the occurrence of dynamical bi¬ 
furcations (Eig. |4]), as it is a non-linear process ( 3 ^ . 

There exist a number of different transport mod¬ 
els that have been developed till now. Among them 
are intranuclear cascade calculations, time-dependent 
Hartee-Fock model, Boltzmann-type kinetic equa¬ 
tions based on hadronic non-equilibrium quantum 
transport field theory, molecular dynamics methods 
(MD, QMD, AMD, FMD), nuclear fluid dynamics 
calculations, etc. The detailed analysis can be found 
elsewhere (e.g., 0100), but, summarizing all pros 
and cons, it is possible to conclude that although 
much progress has been made over the past couple of 
decades, we are still far from having models that are 
formally well founded, practically applicable, and suf- 
hciently realistic to be quantitatively useful. The de¬ 
scription of nuclear fragmentation dynamics requires 


ISSN 2071-0186. Ukr. J. Phys. 2015. Vol. 60, No. 8 


711 












K.V. Cherevko, L.L. Jenkovszky, V.M. Sysoev et al. 


8 


6 

A 

Ll. 

2 

^ 4 
V 


2 


0 

Fig. 5. Mean IMF multiplicity obtained from different statisti¬ 
cal models. The open and solid stars on the right panel are data 
from the central collisions of ^^^Sn + ^^^Sn and ^^^Sn + ^^^Sn 
systems at i? = A = 50 MeV 
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Fig. 6. Predicted mass distributions from five evaporation 
codes for neutron-deficient system 1 (left panel) and neutron- 
rich system 2 (right panel) 14ill 


that the proper account for the basic quantal nature 
of the system be taken 

Similar situation is observed for the later stages 
of the collision. One of the considerable problems is 
that the time needed for the equilibration and the 
transition to the statistical description is still under 
debate. The other difficulty to be overcome is that 
most statistical multifragmentation codes have been 
developed to describe specific sets of data and nearly 
all of them have different assumptions. They are not 


equivalent [4l| . Basically, two main sets of models are 
used to describe the final stages of the reaction that 
are the multifragmentation models and the evapora¬ 
tion models. From the data available in the literature 
(4l| . one can see (Figs. [S] and [5]) that the overall be¬ 
havior of the systems shows a similarity among mod¬ 
els suggesting some predictive power as for system’s 
behavior in the multifragmentation phenomena and 
allowing for suggesting physical processes responsible 
for the formation of fragments. At the same time, the 
behaviors of single observables differ from model to 
model suggesting some limitations on their applica¬ 
bility in treating the experimental results. 

In spite of a long history and a high number of dif¬ 
ferent approaches used, there is still a number of prob¬ 
lems left that have no explanation. For example, the 
models involving a phase transition have the strong 
position based on the recent works dealing with the 
bimodality (^ . But, in the same time, their justifi¬ 
cation is quite questionable. They are insensitive to 
many initial parameters [2l| , when it is known for or¬ 
dinary liquids that the phase transition is usually the 
quantity difficult to be prepared. The most success¬ 
ful models in describing the observed fragment mass 
distribution nowadays are the statistical equilibrium 
models mil that try to reduce the intractable prob¬ 
lem of a time-dependent highly correlated interacting 
many-body fermion system to the much simpler pic¬ 
ture of a system of non-interacting clusters [ 43 . First 
based on Bevalac’s results, they were used for inter¬ 
mediate energies of HIC with many substantial vari¬ 
ations [il. In the same time, such models have a 
number of difficulties 2l|,l43|. The question with the 
equilibrium at a low freeze-out density that lies in 
the basis of such theories is not confirmed by some 
microscopic approaches. There are evidences of only 
a small part of the system being heated illil- The 
other problem concerns with the predicted kinetic en¬ 
ergy of fragments, bein g lo wer than the temperature 
of the gas of fragments 4J|, etc. 

Therefore, the determination of the mechanism 
responsible for the multifragmentation phenomenon 
and able to solve the existing problems in natural 
way is of particular interest as this issue has not been 
settled up to now. 


2.1. Thermodynamics of heavy ion collisions 

The interest to the intermediate-energy collisions is 
triggered by a fact that the multifragmentation is 


712 


ISSN 2071-0186. Ukr. J. Phys. 2015. Vol. 60, No. 8 






















Common Approaches in Description of Ordinary Liquids 



Fig. 7. Left: the pressure-temperature coexistence curve for bulk nuclear matter. Right: the temperature- 
density coexistence curve for bulk nuclear matter. Errors are shown for selected points to give an idea of 
the error on the entire coexistence curve 14611 


claimed to be connected with a phase transition in 
NM (^. Such an assumption comes from the ther¬ 
modynamic analysis of the coexistence curve in NM 
(Fig. [7]) based on the intermediate HIC data. In that 
case, the Fisher droplet model is widely used due 
to the fact that the distribution of fragments can 
be described by Fisher’s well-known formulae for the 
droplet mass distribution F7, 4^ 


dP 

dM 




( 2 ) 


where P is the fragment formation probability, M 
is the mass, and a is the critical exponent. In that 
case, the basic principles of the analysis are absolutely 
the same as in the physics of ordinary liquids. The 
description of the formation of clusters is based [i^ 
on the Gibbs balance equation for the free energy G: 


AG = AE- TAS + PAV, (3) 

where AE and AS are the energy and entropy costs 
of the formation of a cluster, respectively, P is the 
pressure and AV is a change in the volume owing to 
the formation of a cluster. Within such approach, 
it appears next to be possible to use the classical 
Guggenheim relation for the coexistence curve (50| : 


Pl,v , 

-= 1 -I- 0i 

Pc 



±&2 



(4) 


where 6i, &2 are the parameters and /3 = 

is the critical exponent. Equation (|3]) together with 


the experimental vapor branch of the coexistence cur¬ 
ve Ei 


P_ ^ 
Pc 


allow constructing the complete coexistence curve 
of nuclear matter (e.g., see Fig. III). Within this ap¬ 
proach, the modern theory of phase transitions to¬ 
gether with the thermodynamics allow estimating the 
critical temperature in inhnite NM and a finite nu¬ 
cleus. That can give, in turn, a deep insight into the 
properties of NM and can be of great importance for 
developing the proper NM EOS. Quite promising is 
the thermodynamic analysis of the phase trajecto¬ 
ries of systems formed in HIC. Recently, there ap¬ 
peared a work devoted to the study of different possi¬ 
ble phase trajectories of the excited nuclear system at 
the P — V plane [sH. It was aimed to get a qualitative 
picture of the phenomenon considering the boundness 
of the system and the existing laws, which describe 
the behavior of the ordinary liquids in the metastable 
and spinodal states [Hi. Assuming the system before 
the collision to be at point L of the phase diagram 
(Fig. [8]), the different scenarios of system’s evolution 
were studied. 

Evolution of the system crucially depends on the 
excitation energy, mass number, and energy re¬ 
lease conditions. This position suggests two different 
groups of phase trajectories, namely: the single-phase 
transition (marked with the dashed line) and two- 
phase transitions (marked with the solid line). It is 
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Fig. 8. Possible phase trajectories of a nuclear system in the 
proton-induced multifragmentation phenomenon 


obvious that there could be a mixture of two de¬ 
cay channels, when different parts of the system are 
found in different areas of the phase diagram. Thus, 
one may say that the nuclear multifragmentation is a 
nonlinear phenomenon, which means that the quali¬ 
tative picture is different depending on system’s pa¬ 
rameters. The thermodynamic analysis conducted in 
the work has shown that not all of the mechanisms 
could be realized in nuclear systems because of their 
size. Some mechanisms does not meet the require¬ 
ments for the time needed for the process. Summing 
up all pros and cons of different decay channels and 
qualitative characteristics, it was suggested that the 
spinodal decomposition could not be responsible for 
the multifragmentation phenomena because of the 
system size. That fact can be explained from Cahn’s 
theory of separation by the spinodal decomposition 
( 5 ^ . In that theory, changes in the free energy read 


AF = 


1 

2 




-I- K{Wcf 


dV, 


( 6 ) 


where / is the free-energy density of a homogeneous 
material with composition c, K{yc)^ is the addi¬ 
tional free energy density if the material is in a 
gradient in the composition. Therefore, the system 
is unstable relative to the Fourier components with 

( q 2 f j fi ^ 

-sufficiently large wave¬ 
lengths, as this decreases the free energy, when the 
system is in the unstable region. This result shows 
that, for the smaller wavelengths, there is no de¬ 
crease in system’s free energy, and there is no sign 


change in the diffusion coefficients. Therefore, for a 
small enough system, the spinodal behavior cannot 
be justified. As for the most appropriate decay chan¬ 
nel, it appeared to be the mechanical breakdown of 
the shell in a single-phase process that may be fol¬ 
lowed by the nietastable boiling. Although, there also 
could be a mechanism based on the nietastable boiling 
of the inner part of the system. The suggested ther¬ 
modynamic analysis allowed for a simple qualitative 
picture of the phenomena. At the same time, some 
further quantitative estimations based on the com¬ 
bination of macroscopic and microscopic theories, as 
well as on computer simulation results, are needed. 
There exist quite a lot of other examples (e.g.. 


54l | and references therein) of the successful imple¬ 


mentation of statistical thermodynamics approaches 
to quantum systems and nuclear matter, in particu¬ 
lar. This suggests the thermodynamics to be an im¬ 
portant ingredient in the nuclear matter studies. 


2.2. Exotic topologies in the intermediate 
energy collisions. Hydrodynamic approaches 


The hydrodynamic description of nuclear matter 
dates back to the 1980s (see, e.g., [Ellsi) and is 
now widely used for the high energy HIC. In some pi¬ 
oneer works [^, the possibility to use the hydrody¬ 
namic description in the lower-energy limit was con¬ 
firmed, by basing on the analysis of the nucleon mean 
free path that was defined as A = l.d^fm. Later on, 
quite a lot works addressed the problem of a nucleon 
mean free path and the evaluation of the in-medium 
nucleon cross section, with the Pauli blocking being 
considered Still, the results found in some 

more recent publications with the parametrized 
in-medium nucleon-nucleon cross sections from the 
Dirac-Brueckner approach based on the Bonn-A po¬ 
tential give the values A ^ 1.4 fm and A ^ 1.3 fm 
for p/pq = 1 and p/po = 1-5, respectively, at A = 
= 50 MeV/nucleon. Those results confirm well the ar¬ 
gumentation in (^ . Despite that, the hydrodynamic 
description of the HIC at intermediate and low ener¬ 
gies has been the subject of only few theoretical stud¬ 
ies [^. During the last decades in a number of works, 
the possibility of the formation of toroidal (Fig. [9]) 
and bubble (Fig. ITOll structures in the head-on HIC 
was studied extensively with the help of microscopic 
transport models . Within the Boltzmann- 

Uehling-Uhlenbeck (BUU) model, it was shown that. 
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at the energy interval 60-75 MeV/A, the different 
exotic topologies of the formed objects can be ob¬ 
served, depending on the stiffness of the nuclear 
mater EOS. Main peculiarities observed in the qual¬ 
itative picture, which comes from the Boltzmann- 
Uehling-Uhlenbeck calculations and the available ex¬ 
periments [silil, are the fact that the expanding 
velocity of the outer surface of the system is much 
smaller than the expanding velocity of the inner sur¬ 
face, simultaneous breakup with few fragments of 
similar masses and low kinetic energies and with 
the angular distribution of fragments that is almost 
isotropic in the case of soft EOS (the incompressibil¬ 
ity at saturation density Kq = 200 MeV), and the 
presence of fragments in the plane perpendicular to 
the beam direction for the stiff EOS (the incompress¬ 
ibility at the saturation density Kq = 380 MeV). For 
some time, the experimental results were not able to 
confirm the occurrence of the predicted geometries 
[g^. But, later on, the signatures of the “doughnut”- 
like structures with the production of similar-size 
intermediate-mass fragments were observed in cen¬ 
tral HIC and again confirmed by transport mod¬ 
els 0 . Unfortunately the calculations with trans¬ 
port theories used for the phenomena in focus are 
not able to give reliable information on multiparticle 
observables at the late stage of the process due to 
the fact that they do not include multiparticle corre¬ 
lations and fluctuations [^. In some works ^ 671, 
the Rayleigh-Taylor mechanism was suggested to be 
responsible for the formation of fragments, when the 
exotic topologies were explained with the help of 
shock waves. Such type of structures is also some¬ 
what similar to those studied theoretically from an¬ 
other approach [^. Different shapes observed in the 
experiment may yield the valuable information about 
the nuclear incompressibility, which is a key param¬ 
eter of the nuclear EOS, as well as the informa¬ 
tion regarding the surface tension and NM viscos¬ 
ity. The possibility for different shock wave mecha¬ 
nisms to be involved into the process can also give 
an extra link in between the Boltzmann-like trans¬ 
port theory and the hydrodynamic approach. Re¬ 
cently, it was shown [69| that the observed quali¬ 
tative picture has much in common with the col¬ 
lisions of high-speed ordinary droplets 73-72|. For 
the latter in the certain energy range, the forma¬ 
tion of “doughnut”-like structures with the follow¬ 
ing fragmentation into several secondary drops of ap- 
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proximately equal masses is observed l73ll. w hich has 
much in common with the nuclear case |64| . In a re¬ 
cent work [g^ . the system of two identical heavy nu¬ 
clei (e.g. ®^Nb + ®^Nb) involved in a head-on col¬ 
lision was considered. The symmetry of the system 
allowed the simplification of the model by chang¬ 
ing to the collision of a spherical nuclei of radius 
R and density po with a rigid wall that moved 
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Fig. 11. Different stages of the evolution: symmetry plane im¬ 
pacting the nuclei a; lateral jetting, when the shockwave veloc¬ 
ity becomes equal to the contact edge velocity b. Compressed 
zone is shown in grey. The geometry of the system in the case 
of “stiff’ EOS {K = 380 MeV) c; geometry in the case of “soft” 
EOS {K = 200 MeV) d 


toward it with the velocity vq (Fig. [TTl a). The 
slip boundary condition was applied on the wall 
surface to account for the difference in a viscous 
behavior. 

From the analysis of the recent theoretical and ex- 
perimental results for the collisions of liquid droplets 
together with the BUU 641 and nuclear fluid 
dynamics (NFD) results for HIC, it is suggested 
that there should be four distinct stages of the colli¬ 
sion [Fig. [TT]. Among them are the violent stage at 
the beginning of the process, when the highly com¬ 
pressed zone is formed, and the second stage that is 
characterized by the lateral jetting and lasts till the 
shock wave from collision plane reaches the boundary 
of the nuclei. During those two stages, the final topol¬ 
ogy of the system is fully defined. At the later times, 
the system goes through the third stage correspond¬ 
ing to the expansion process and finally comes to the 
last stage, when the fragmentation takes place. The 
main quantitative estimations within the suggested 


Quantitative characteristics 

of the process for the system 

at 60 MeV/nucleon with “stiff’ EOS {K = 380 MeV) 



m 

[63,^ 


Reaction time ra, fm/c 

118 

120-160 


Maximum density p/po 

1.6 

~1.5 


Spreading size — 

2.6 



Rim radius at ra, fm 

1.9 



Number of similar mass IMF 

5.6 

3-6 

4-6 


model and their comparison with the BUU calcula¬ 
tions and experimental data are presented in Table. 

One can see that the above model slightly over¬ 
estimates the spreading size and underestimates the 
reaction time. This probably comes from the ap¬ 
proximations used in the model. For example, ac¬ 
counting for the viscous dissipation and changing 
to the Navier-Stokes equations, as well as account¬ 
ing for the changes in the surface tension coeffi¬ 
cient due to the temperature increase, can influ¬ 
ence both of the above values. The model doesn’t 
give the strict results, but rather provides a picture 
of the involved physical processes as clear as pos¬ 
sible. For such a model with no adjustable parame¬ 
ters, the results for the maximum density and overall 
timescale of the process are in good correspondence 
with the BUU calculations, ft is also worth mention¬ 
ing that the result for the number of fragments of 
similar masses is in good correspondence both with 
the experiment and the BUU calculations. Such re¬ 
sults show that the hydrodynamic description based 
on the laws developed for ordinary liquids seems to 
explain the overall behavior of the nuclear systems 
in the head-on HIC in focus and allows a simple 
physical picture of the formation of exotic struc¬ 
tures. 

One can conclude that the combination of the hy¬ 
drodynamic approach together with the transport 
theory calculations can reveal the physical nature of 
the phenomena observed in HIC and give possibility 
to proceed with extracting the data on the NM prop¬ 
erties from the HIC at different energies. 

3. Surface Effects in Nuclear 
Matter: Macroscopic Description 

The nuclear matter properties have been extensively 
studied for more than half a century. There exist quite 
a number of theoretical works devoted to different 
models of inhnite nuclear matter, as well as hnite nu¬ 
clei. When discussing the different approaches, one 
can divide them into two groups that are microscopic 
and macroscopic approaches. At first glance, the mi¬ 
croscopic approaches seem to be much more general 
than the macroscopic ones that deal only with the 
averaged values and can be obtained basically as the 
average of the microscopic theory. At the same time, 
the exact microscopic calculations are very challeng¬ 
ing to perform, and the precise data on the micro- 
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scopic parameters are needed that makes it difficult 
for practical applications. During recent decades, the 
impressive progress has been achieved in the macro¬ 
scopic description of nuclear matter (d^ . [zllzi- A 


number of papers devoted to the thermodynamics 
of small systems or the hydrodynamics of nuclear 
matter appeared [i, |5l[53,[zi,[83- Therefore, the nu¬ 
clear systems that consist of a number of nucleons 
that is neither very small nor very large seem to be 
a good object for being studied with the help of a 
combination of both macroscopic and microscopic ap¬ 
proaches. The bright example of such a cooperation 
of the methods in nuclear physics is the liquid droplet 
model [ 8 l| introduced in the 1960s, which makes pos¬ 
sible the description of the average properties of a 
saturated system such as a nucleus consisting of two 
components (neutrons and protons) with account for 
the boundary effects and the presence of a diffuse 
layer. In that case, the energy at equilibrium is given 
as 

E = (^-ai -h J6^ + A + 

+ 0,2(1 + 2 e)A.^ + Qt ^+ a^A^ H- 
Z‘^ f _ 1 ,_A 


— C2Z^A3 — C3—;-C4 


Z3 


(7) 


with 


= 1 - 3e, 


P_ 
Po 


( 8 ) 


where J, K, M, Q, ai, 02 , and 03 are the coefficients 
evaluated from the observed averaged phenomenolog¬ 
ical inputs and some relations based on the differ¬ 
ent microscopic considerations, A is the nuclear mass 
number, and Z is the charge number. As a result, 
the liquid droplet model is able to provide a quite 
precise description of the properties of finite nuclei. 
As comes from the liquid droplet model in studying 
the curvature-correction term for the nuclear matter, 
one should keep in mind the connection between the 
surface and bulk properties of the matter 81, 8^ . As 
shown in the droplet model, the coefficients in the 
term proportional to As in the expansion of nuclear 


properties in terms of the fundamental dimensionless 
ratio A~3 are related to the bulk properties of NM 
described by the terms proportional to A and As. 

Corrections due to a curvature may play an im¬ 
portant role, when studying light nuclei or processes, 
where surface terms are important. Particularly im¬ 
portant are those corrections in the interpretation of 
multifragmentation experiments 8 ^ 8 ^ - 8 ^ , in which 


fight nuclei necessarily appear. The exponential de¬ 
pendence of the yield of fragments on the surface 
tension makes this process sensitive to the curvature 
corrections 4^8^. Other important phenomena that 
may be affected by changes in the surface tension due 
to curvature corrections are: 

a) the appearance of the neck region in the fission 
processes and 

b) the hydrodynamic instability of the structures 
formed in heavy ion experiments governed by surface 
effects (69l . l8^ . 

Quite a number of papers devoted to studies of the 
surface energy and the properties of surfaces in NM 
iiSi lave appeared. Furthermore, the dependence 
of the surface tension coefficient (and surface energy) 
on the surface curvature and its influence on differ¬ 
ent physical properties were also studied by various 
groups of researchers 8^l91|. Still for decades, it re¬ 
mains a controversial issue in mesoscopic thermody¬ 
namics (q^ - I^ . 

The thermodynamic description of the curvature 
correction, originating from the difference between 
the equimolar surface and the surface of tension 
91 ,9^, dates back to the 1940s. The Tolman length S 
was originally introduced in to describe the cur¬ 
vature dependence of the surface tension of a small 
liquid droplet. It was defined as a correction term in 
the surface tension a of the liquid-vapour droplet in 
the isothermal case: 


r(i?) = 


1 - 


R 


(9) 


where R is the droplet radius equal to the radius of 
the surface tension |95l . l96j . and CToo is the surface ten¬ 
sion of a planar interface. Eq. (O originates from the 
Gibbs-Tolman-Koenig-Buff’s thermodynamic equa¬ 
tion and the assumption that 6 is independent of 
R for 5 <C i? [^. This physics should work not 
only for liquid droplets, but also for any system with 
curved interface of a iioii-negligible boundary layer 
[ 9 ^. This situation corresponds to nuclei and nuclear 
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Infinite nuclear matter 
o ° o o o o 


Nucleons external to the 
chosen volume are removed 



matter 


ject for continuous debates with implications in the 
early Universe and relativistic HIC. The interest in 
this subject is triggered by the interpretation of the 
observed elliptic flow in heavy-ion collisions at RHIC 
Tni - [Tnfil | as a manifestation of the low shear viscos- 


ity o f nuclear matter, and a subsequent realization 


Imaginary surface of the ^ Real surface of the 
chosen volume with A nucleons shrunk 

Fig. 12. Schematic picture of the gedanken experiment l99ll 

systems with a finite diffuse layer Iz3- Based on that 
fact, a thermodynamic approach that links the EOS 
of NM with the curvature corrections to the sur¬ 
face tension coefficient has been introduced recently 
[9^. In this work, the analysis is based on the simple 
gedanken experiment. First, one should consider infi¬ 
nite nuclear matter [Pq, T = const) with the chosen 
spherical volume Vq = in it consisting of A 

nucleons. If all the nucleons outside the chosen vol¬ 
ume are removed, one gets a “nuclear droplet” that, 
due to the surface tension, shrinks to the volume 
V = where R is the final radius of the cho¬ 

sen volume (Fig. 1121) . This droplet remains in equi¬ 
librium. The analysis of changes in the surface size 
due to the shrinkage allows calculating the curvature 
correction to the surface tension coefficient from the 
EOS of NM. The Tolman length S is calculated for 
different effective interactions (SLy6, SkM* and SV- 
min). The calculated results for S in the case of, for 
example, the S'U-min parametrization, appeared to 
be 5 ^ —0.55 fm that is close to the distance be¬ 
tween the nucleons. That result correspond to the 
theoretical estimations based on the Gibbs-Tolman 
formalism. Therefore, the statistical thermodynamic 
approach can be useful in studies of the surface prop¬ 
erties of nuclei. 

4. Hadronic Matter and Quantnm Liquids 

Nowadays, there exist quite a lot of papers devoted 
to the application of classical results of the phase 
transition theory and the theory of ordinary liquids 
for hadronic matter and quantum liquids. That can 
be explained by the growing interest in and the un¬ 
derstanding of the universality of critical phenom¬ 
ena in different media over a tremendous span of 
substan ces - fro m ultracold atoms to dense nuclear 

ions of the nucleation rate 


108 


1071 of t he universality of this phenomenon, see also 
Most remarkably, these studies penetrated 


109 


110l - lll4j different neighboring fields of physics such 


as condensed matter and statistical mechanics and 
their interface with nuclear and astro-particle physics. 

4.I. Quantum liquids 

When applying the classical fluid thermodynamics 
to nuclear matter, quantum corrections become im¬ 
portant as compared to the motion of atoms and 
molecules in ordinary liquids 115 . E3. There is 
an essential difference between nucleons, on the one 
hand, and nuclei, atoms, or molecules, on the other 
hand. In the latter case, the relevant energies - bind¬ 
ing energies, exciting energies, etc. - are much smaller 
than the rest frame masses of the system or of the con¬ 
stituents. In the case of hadrons, they are of the same 
order as the mass of the system. Hence, the quantum 
effects are of importance. 

The amount of quantum corrections to the classi¬ 
cal theory can be quantified from the expansion in 
powers of the Planck constant h of the quantum sta¬ 
tistical sum of a system of interacting particles. For 
the H elmholtz free energy, the following expression 
holds ins Ill7l |. in which the first quantum correc¬ 
tion is proportional to : 


R — Rclass 


1 + (A) ujiTwn 


( 10 ) 


Here, Fliass is the classical free energy, oj is some di¬ 
mensionless function of the order of 1, and a new 
dimensionless parameter A is defined as 


A = 


aWme 


( 11 ) 


during the (de)confinement phase transition are sub- 


where m is the mass of the particle. A is called the de 
Boer parameter and describe s the amo unt of quantum 
effects in statistical systems fnR fIl7l | . 

The apparent similarity of the properties of differ¬ 
ent media sug gest s the use of the law of corresponding 
states (LCS) |ll5l | . For real atoms and molecules, the 
intermolecular potential $(r) can be reduced, with 
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reasonable precision, to a product of the energy con¬ 
stant e (the depth of the minimum in the particle 
interaction potential) and a function of the dimen¬ 
sionless variable x = 

<^{r) = e(l){x), (12) 

where a is the position of the minimum. Then, by 
LCS, for a group of substances with a common po¬ 
tential 4>ix), the universal functions for the pressure 
p* = p* energy E* = E* (T*, V*), free en¬ 

ergy F* = F*{T*,V*), etc. should exist, where 
T* = — is the reduced temperature and V* = ^ 
is the reduced volume. As a consequence, the ther¬ 
modynamic properties of all the substances from a 
given group are similar in terms of a properly chosen 
scale. In view of the similarity of the nuclear forces 
to the van der Waals forces, see Fig. [TI] LCS should 
work both in ordinary liquids and in nuclear matter. 

Recently, it was suggested to extend the principle of 
corresponding states for ordinary liquids to the pos¬ 
sible fluid-superfluid and conductor-superconductor 
transitions for a variety of systems, ranging from 
^He, "‘He, and superconducting electron fluids to nu¬ 
clear matter. Such an extension is based on the de¬ 
scription of the superfluidity phenomenon in terms of 
wave packets. Superfluidity assumes the appearance 
of phased waves (wave packets) of bosonic quantum 
particles due to the Bose condensation (Fig. [13]). 

As a consequence, it appeared to be possible to 
define the analogue of the d e Bo er parameter - di¬ 
mensionless number BJEMS [llSl ] 

BJEMS = g I /A (13) 

Analogously to the de Boer parameter, it takes the 
quantum effects into account. The difference is that 
it considers the exchange effects. Therefore, it might 
be used for the superfluid phase transitions. The pa¬ 
rameter BJEMS is assumed to be of the same order 
of magnitude for all fluid-superfluid and conductor- 
superconductor transitions. This hypothesis is based 
on the general principles, rather than on the micro¬ 
scopic theory, and it appears to be correct for all the 
real Fermi and Bose systems from the same univer¬ 
sality class. 

J 1 .. 2 . Hadronic matter 

Studies of the hadronic matter are an important part 
of the modern astrophysics and HIC physics. In the 



Fig. 13. Distribution of waves of particles: in a classical sys¬ 
tem (high temperatures) (a); exactly at the transition temper¬ 
ature (6); slightly below the transition point (c) 


early Universe, a particular nucleation dynamics dur¬ 
ing the confinement can influence the inhomogeneities 
in the hadron distribution. Nowadays, the possibility 
for a number of colored objects such as quarks and 
gluons to escape the hadronization du ring the c ooling 
of the Universe is intensively studied [mlmi . This 
requires the understanding of the surface properties 
of the formed quark “nuggets”. Statistical thermody¬ 
namics can be quite helpful in the studies of those 
properties. In HIC, superheated hadron matter can 
produce extra entropy and might lead to observable 
consequences, detectable in the experiment. Earlier, 
there were the attempts to use statistical thermody¬ 
namic approaches for studying the hadron matter, 
confinement, phase transition s, and early evolution 
of the Universe 8^ 122l4l24l |. Within the suggested 
approach, some interesting results were obtained in 
the description of the expansion of the Universe dur¬ 
ing the early stages of its evolution. The behavior of 
the surface tension coefficient in the (de)confinement 
phase transitions was studied. The obtained results 
for a parton gas can be useful in the studies of the 
color-glass condensate that is discussed in the frame¬ 
work of the studies of relativistic HIC. 

The similarity of the behavior of ordinary liq¬ 
uids and hadronic matter, as well as the exis¬ 
tence of a number of interesting results obtained for 
hadronic matter from the statistical thermodynamic 
approaches, suggests the application of classical laws 
in high-energy physics. 


5. Summary 

Macroscopic theories well-developed for ordinary liq¬ 
uids have shown themselves to be successful in 
describing the average properties of nuclear sys- 
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terns. They work well for nuclear matter in combina¬ 
tion with microscopic approaches naturally comple¬ 
menting each other. Hydrodynamic approaches de¬ 
veloped for ordinary liquids are widely used in the 
studies of heavy ion collisions at different energies, 
as they explain well the observed phenomena, and 
their results can be easily compared to the calcula¬ 
tions within transport theories. The thermodynamic 
description has shown to be fruitful for finite nuclei 
and intermediate-energy heavy ion collisions. 

At the same time, statistical thermodynamics can 
be applied not only for ordinary liquids and nuclear 
matter, but also for studies of the properties of quark- 
gluon plasma. It has shown to be useful in under¬ 
standing the behavior of quantum liquids in phase 
transitions. 

Macroscopic theories can be used in constructing 
a proper equation of state of nuclear matter that is 
one of the long standing and desirable goals in nu¬ 
clear physics. They can be used in the studies of the 
intermediate heavy-ion collisions, which are a window 
into the nuclear matter equation of state and can be 
treated as the analogue of the P-V-T studies in the 
physics of ordinary liquids. 

Statistical thermodynamics and phase transition 
theory widely used in the physics of ordinary liq¬ 
uids work well for hadronic matter. They are quite 
fruitful in studying the (de)confinement phase tran¬ 
sitions. Interesting results are found in the works de¬ 
voted to the color glass condensate that is discussed in 
the context of relativistic heavy ion collisions. They 
can be useful in studies of the early evolution of the 
Universe. 

As a summary, one can conclude that the ap¬ 
proaches originally developed for ordinary liquids are 
widely used in the studies of hadronic matter and can 
give valuable results. 
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K.B. HepcGKO, JI.JI. CHKoecbKuu, 

B.M. Cucoee, OeH-IlIoy ?Kam 

SAFAjiBHHH ni;i,xi;i, 

AJJH OHHCy KJlACHHHHX PIAHHHHX 
CHCTEM 1 HAEPHOI MATEPH 
P e 3 lo M e 

B poboxi pobHTbcn cnpoba KopoTKoro orjin^ MOi«;jiHBOCTeH 
SaCTOCyBaHHH CTaXHCTHUHOi’ XepMO;],HHaMiKH npH BHBUeHHi 
a;],poHHoi" Maxepii’. Po3rjm;],aexBcn: Mo:«;jiHBicxB 3acxocyBaHHH 
ri;],po;i;HHaMiKH aJih onncy 4)OpMyBaHHn pi3HHx xonojioriuHHx 
cxpyKxyp B jioboBHx 3ixKHeHHn:x Ba:>KKHx ioniB 3a Majinx enep- 
rin. HoKasano, u];o xaKHH ni/T,xi/T, ^03BOjme oxpnMaxn anajii- 
XHHHHH onnc KOiKHoi" cx a.ni V npoii;ecy xa BHanauHXH BJiacxHBO- 
cxi n:;],epHoi‘ Maxepii" (noBepxneBHH naxnr, cxHCXjiHBicxb xa in.) 
3 aHajii3y BjiacxHBOCxeii yxBOpennx 4>parMeHxiB. Po3rjm;],aio- 
xbcn nepeBarn xepMO^HnaMiuHoro anajioy (J)a30BHx xpaexxo- 
pin cncxeM, m,o yxBopioioxbcn b 3ixKHeHHHx Ba^cKHX ioniB. Po3- 
rjin^aexbcn MoncjiHBicxb BHKOpHCxanHn cxaxncxHnnoi xepMO- 
;i;HHaMiKH AJiH BHanaueHnn nonpaBOK ;],o Koe4)iLi,ieHxy noBepx- 
neBoro naxnry, BHKJiHKannx kphbh3hoio noBepxni, 3 piBnnHnn 
cxany n^epnoi Maxepii. Po3rjmHyxo npnKJia^H 3acxocyBaHHH 
cxaxHCXHHHoi xepMo;i,HHaMiKH B ^ocjii,ip«;eHHnx a^ponnoi Ma¬ 
xepii xa KBanxoBHx pi^nn. 
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